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HEAT  TRANSFER  TO  A  SPHERE  AT  THE  TRANSITIO^.  FROM 

FREE  MOLECULE  FLOW 

P.  Hammerling  and  B,  Kivei 


SUMMARY 

Baker  ana  Charwat^  have  recently  published  a  semi-quantitative 
study  of  the  correction  to  the  drag  of  a  sphere  at  the  breakdown  of  the 
free  molecule  flow  region.  Their  expression  for  the  drag  coefficient  it:< 


which  shows  the  dependence  on  flight  velocity.  V^,  mean  emitted  particle 
velocity,  V^,  body  radius,  a,  ambier.t  number  density,  p^,  and  collision 
ci  oss-section,  S 

eo 

Using  Similar  assumptions  to  those  in  the  above  work  but  slightly 
different  methods,  we  find  the  total  energy  given  to  a  sphere  is 


where  is  the  ambient  mass  density.  Our  coefficient,  0.2,  is  smaller 

than  the  corresponding  coefficient  in  the  drag  expression  0.4c.,  found  by 

Baker  and  Charwat,  The  coefficient  t*.  2  is  estimated  in  a  manner  described 

below,  lii  place  of  the  momentum  recoil  term,  0.  44  V  /V  ,  there  is  the 

e'  o 

factor  i  -  fo  account  for  the  ene.-gy  carried  away  by  the  emitted 

particles.  The  correction  to  the  scattering  distribution  for  the  velocity  of  the 


‘R  m  l.  Baker,  Jr.  and  A.  F  C.harwat,  Phys.  of  Fluids,  1,  73  (1958). 


erni^teo  par*^icie»  not  being  negligible  compared  to  the  flight  speed  has  net 
been  made. 

Our  rTiOdei  starts  from  the  observation  that  the  angular  distribu¬ 
tion  of  emitted  particles  has  simple  analytic  forms  in  the  extremes  of 
being  very  close  to  and  v>'>ry  far  from  the  sphere.  At  the  surface  the  density 
of  emitted  particles  lor  diffuse  scattering  is  2  co“8/V^.  The  angle  6 

IS  made  by  the  .'me  of  flight  and  the  line  connecting  the  center  of  the  sphere 
ard  the  point  of  interest.  Tlie  factor  cos  6  occurs  because  of  the  decrease 
at  6  of  the  projected  area  normal  to  the  flight  direction  per  unit  surface 
area  of  the  sphere  The  factor  2  arises  because  of  the  Lambert's  law 
scattering  Far  from  the  surface  the  distribution  can  be  obtained  by  count¬ 
ing  all  particles  emitted  with  the  same  direction.  At  the  distance  r  this  is 


oC’ 


Co< 


c^ 


_! 


(3) 


in  agreement  with  the  expression  given  in  the  paper  by  Baker  and  Charwat. 

Another  difference  lies  in  the  approxima  'on  used  to  treat  the  probability 

of  scattered  particles  returning  to  the  sphere.  Near  the  sphere  we  use 

2  1^1  -  {l/2i  sin  u  (a/r'-  per  collision,  and  at  larger  radial  distances. 

? 

2  cos6  (a/r)^  The  .mrmer  is  chosen  h*'  :ause  for  a/r  =  1  at  6  =  0,  both 
scattered  particle.-  hit  the  sphere  while  *t  8  =w/2  just  one  of  the  scattered 
particles  hits  the  sphere  However,  the  result  does  not  seem  sensitive  to 


this  difference 

We  have  used  the  cos  t  angular  distribution  with  a  (a/r)^  dependence 
out  to  t-AO  radii  from  the  sphere  surface  At  larger  radial  distance^ 

Eq  /3‘*  is  used  With  these  distributions  we  are  able  to  estimate  the  net 
loss  of  energy  by  scattering  It  is  convenient  to  compare  the  total  energy 


brought  to  the  sphere  by  scattered  particles  (q  .  )  with  the  energy  incident 

ga  xH 

particles  that  are  hit  wc-old  have  given  to  the  sphere  (q,  ).  Since  some 

xO  s  s 

of  the  part5>}.es  incident  are  deflected  only  slightly  so  that  they  still  hit  the 

sphere,  it  is  necessary  to  point  out  that  such  events  are  counted  twice. 

They  are  included  as  a  total  loss  in  their  contribution  to  the  sphere 

is  added  to  q  .  .  This  is  used  instead  of  giving  the  net  loss  to  q,  and 
^gain  ^  ®  noss 

no  ad,  jstment  to  q  .  .  We  find  that  the  ratio  q  .  /q,  s  0.  85  from 
'  gain  ^gain  ^ioss 

both  the  near  and  the  distant  regions.  Because  this  ratio  is  the  same  for 
both  regions  our  extrapolation  to  intermediate  distances  while  crude  may  be 
called  semi-quantitative.  It  is  in  this  way  that  we  estimate  the  factor  0.  2. 

The  corresponding  coefficients  0.  2  and  0.  46  should  in  principle  be 
the  same.  In  order  to  compute  the  mcmentam  component  in  the  flight 
direction  and  the  energy  that  can  be  transferred  to  the  sphere  by  a  scattered 
particle  one  uses  the  incident  values  of  these  quantities  multiplied  by  the 
same  quantity,  cos^Q  ,  Admittedly  both  models  are  semi-quantitative. 

That  of  Baker  and  Charwat  does  not  emphasize  the  region  close  to  and  in 
the  limb  of  the  sphere.  Our  own  extrapolation  to  intermediate  distances 
is  open  to  criticism. 

If  we  calculate  the  number  of  particles  hitting  the  sphere,  we  find 
a  net  gam,  whereas  there  is  a  net  energy  loss.  This  results  from  the 

cos^h  term  in  the  energy  return  expression,  which  greatly  reduces  the 
amount  of  energy  return  near  6  =  v/Z. 
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Unfortunately,  the  uni,. ^  '•tainty  m  the  coefficient  is  magnified  as  it  results 
from  the  small  net  betwei.r.  loss  and  gain  events.  It  would  be  worthwhile, 
therefore,  to  see  this  problefn  solved  with  more  accuracy. 

It  is  interesting  to  note  the  important  role  of  surface  temperature 
in  these  results  Jf’or  lower  surface  temperature  the  emitted  particles 
move  away  more  slowly.  By  continuity  their  density  is  increased.  Also, 
their  mean-free  path  is  decreased  This  makes  scattering  events  occur 
closer  to  the  surface  and  leads  to  greater  breakdown  of  free  molecule  flow. 

At  100  miles  altitude  a  20  cm  radius  spherical  satellite  with 
velocity  26,000  ft. /sec.  and  with  surface  temperature  300°K  experiences 
1%  reduction  in  the  total  heat  received.  It  should  be  noted  that  this  is 
a  net  effect  and  that  for  a  1%  reduction  in  heat  transfer  about  6%  of  the 
particles  incident  to  the  sphere  are  actually  scattered  by  collisions.  The 
loss  IS  considerably  less,  since  many  of  the  deflected  particles  still  hit 
the  sphere  and  many  that  would  not  have  hit  are  deflected  into  the  sphere. 

We  are  grateful  to  R.  F.  Probstein  of  Brown  University, 

A  Kantrowitz,  and  N  Kemp  of  AVCO  Research  Laboratory  for  sugges¬ 
tions  and  discussions - 

Or.  the  ensuing  pages  we  give  the  details  of  our  treatment  which  have 
led  us  to  the  above  conclusions 
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A  Model  for  a  Ranfied  Gas  Flow  About  a  Sphere 
^  P.  Hammerling 

LNTRODUCTION  AND  AigS&HPTIONS 

In  free  molecule  flow  t**  effects  of  interiViOiecular  collisions  can 
be  ignored  compared  with  collisions  uith  a  wall  or  an  obstacle.  This 
implies  that  the  mean-Jree  path  in  the  gas  much  greater  than  a  charac¬ 
teristic  length  of  the  booy.  We  shall  consider  the  high  speed  flow  of  a 
Tariffed  gas  about  a  sphere  In  such  a  flow  the  reflected  molecules  are 
moving  relatively  slowly  and  form  a  region  near  the  body  whose  density 
is  greater  than  the  free  stream  density.  In  the  usual  treatments  of  free 
molecule  flow  problems^  one  assumes  that  the  reflected  particles  do  not 
change  the  distribution  function  of  th«  incident  gas.  However,  the  presence 
of  the  denser  layer  near  the  body  introduces  the  possibility  of  collisions 
with  the  incident  stream  which  can  alter  the  expected  momentum  and 
energy  transfer.  As  a  first  step  one  might  look  for  the  consequences  of 
specifying  that  the  incident  and  reflected  particles  make  only  one  collision 
with  one  another.  This  restricted  problem  has  been  solved  for  the  case  of 
a  flat  plate  by  Heineman^  and  by  Lunc  and  Lubonski^. 

The  case  of  a  sphere  presents  difficulties  for  two  reasons:  (IJ  The 
attenuation  of  the  reflected  particles  due  to  collisions  introduces  a  compli¬ 
cated  exponential  factor  in  the  integ.ations,  (2)  The  geometry  involved 
also  makes  for  difficult  integrations.  In  order  to  gain  some  understaiid- 
ir.f  of  this  problem  without  first  having  to  do  ail  the  work,  we  have 
attempted  to  simplify  it  even  further.  The  guiding  spirit  of  the  calcula¬ 
tion  is  for  it  to  be  simple  and  yet  include  the  physics  of  the  problem.  We 


Heineman,  Commun  Apl,  Math,  259  (1948) 

Lunc  and  J  Lubonski.  Archiwam  Mechanki  Stoswanej.  4,  597  (1956) 


-5- 


h  -e  spilt  things  into  two  regions,  the  near  and  far  zones.  Far  away  from 
t.^e  sphere  the  geometrical  difficulties  can  be  overcome;  close  to  the  sphere 
»^he  geometry  must  be  considered  more  carefully,  but  the  attenuation  can 
be  igro^ed.  Since  the  results  in  these  two  regions  should  go  continuously 
into  *i;e  another,  we  have  made  an  estimate  of  the  in  between  region  which 
is  protab! y  good  to  within  20%. 

Behind  the  calculations  are  certain  assumptions  which  we  now  list: 

(a)  Collisions  are  elastic  and  are' visualized  as  between  rigid 
spheres,  the  scattering  being  isotropic  in  the  center  of  mass  system  of 
the  two  particles. 

(b)  Particles  which  hit  the  surface  are  emitted  in  a  random 
direction  with  a  probability  proportional  to  the  cosine  of  the  angle  between 
the  normal  to  the  surface  and  the  direction  of  emission.  This  is  known 

as  diffuse  reflection  One  can  regard  the  emitted  particles  as  an  effusive 

stream  coming  from  an  equilibrium  gas  inside  the  sphere. 

p  _  0 

(Cf  The  accomodation  coefficient  a  =  “■ - is  taken  as 

P. 

'-It 

unity  Ej  =  incident  energy  flux,  Ej.  =  reflected  energy  flux,  =  reflected 
energy  flux  corresponding  to  the  body  temperature. 

id)  The  thermil  velocities  of  both  the  emitted  particles  and  the 
gas  proper  are  small  compared  with  the  flow  velocity, 

RESULTS 
A  Far  Pone 

Cc!^' aider  the  sphere  m  the  coordinate  system  shown.  The  2  axis 
IS  ir  t:  •,  direction  of  moi..-n  of  the  sphere  With  respect  to  a  fixed 
dsrc'>.*^ion  ^  in  the  X  2>  plane  a  point  on  the  sphere  will  have  a  polar  angle 


a.  An  observer  at  some  distance  from  the  sphere  will  see  the  reflected 
particles  as  coming  in  a  parallel  beam  from  that  part  of  the  sphere  visible 
to  him.  In  what  follows  we  shall  have  to  integrate  over  the  surface  of 
the  sphere  seen  in  this  way.  If  one  uses  spherical  polar  coordinates  about 


2  the  surface  integrals  will  be  of  the  form 


2 


The  total  outward  flux  in  a  given  direction  0  coming  from  that 
part  of  the  sphere  that  is  visible  can  now  be  found.  This  number  is  equal 
to  the  number  coming  in  at  a  given  point  multiplied  by  the  probability  of 
emission  into  the  0  direction  and  integrated  over  the  visible  surface. 

d  T"  (  0  !  ^  M  U  a 

nr  r  (5) 


/ 


w 


t( 
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The  origin  of  cos  p  in  equation  (6)  may 
be  seen  in  Fig.  {?.).  Consider  the  sphere 
and  a  concentric  shell.  The  velocity 
vector  of  the  particle  emitted  in  the 
direction  a  makes  an  angle  p  with  re¬ 
spect  to  the  normal  to  the  shell.  In 
computing  fluxes  we  need  cos  p  as  in 
equation  (5),  since  we  are  interested 
in  particles  crossing  normal  to  the 

2  2  1/2 

surface.  Using  the  law  of  sines  we  have  cos  p  =  (1  -  (R/r)  Sin  a)  ' 

For  large  r,  cos  P  =  1  which  is  the  case  under  discussion.  The  far  zone 
results  should  be  good  from  a  few  sphere  radii  on  out. 

Since  we  require  the  emitted  particles  to  go  in  the  diiection  Q 
the  emission  law  is  cos  a/ir.  In  our  coordinate  system 


cK.  s  r  :rzos  j  S iniQ  ^  cos  S'cgsQ 


(6) 


The  out^ward  flux  is  then 


9 

z'  ~r!T_  IT 

I  j  ^  CcS  oc  J  S' -h  f  f Sm  S-Cosfi 

where  cm  6  =  tan  cos  6  and  cos  a  is  given  by  {6I1.  This  turns  out 

to  be 


f  (&)  ^  NuTrr'— -  TT  (  q) 

^  3  7r-^  / 

ZL  (  0)  =  (^  ~  =3  )  -  O  -f-  ?  'n  0 


(8) 
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The  total  outward  flux  is  found  by  integrating  12(0)  over  all  6  and 
^  Since  some  particles  are  emitted  in  directions  S’  7  Tl  ;  this  results 


N  u 


The  emitted  particle  density  at  a  distance  r,  Nj^(r)  is  obtained  by 
noting  that  flux  =  Njj(r)Ut,r  d-Zl  , 

hence  flux/soiid  angle  ~  7** (  Q)  =  so  that 


'  ^  "  rr  'XTT^  r=*- 


Due  to  the  possibility  of  collision,  a  pai  '  cle  leaving  the  surface 
will  only  have  a  certain  chance  of  surviving  unhit  to  a  distance  r;  this 
probability  is  given  by: 


P=  e 


,  where  is  the  mean-free  path.  Lj^ 
^  ~  *  when  R- 


With  these  preliminary  ideas  in  mind 
let  us  return  to  a  physical  consideration 
of  the  problem  at  hand.  The  incident 
particles  may  collide  with  the  reflected 
particles  in  such  a  way  as  to  send  one 
of  them  into  the  sphere.  The  proba¬ 
bility  for  such  an  event  must  be  calcu¬ 
lated.  On  the  other  hand,  some  of  the 
incident  particles  that  could  have  hit 
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amcimt  of 


the  sphere  will  not  be  able  to  do  so  because  oC  collisions.  The 

energy  per  particle  that  can  be  carried  to  the  sphere  remains  to  be 

calculated  also* 


Consider  the  cross-section  at  the  left. 
W*‘  shall  consider  all  collisions  within 
the  cylinder  as  losses,  whereas  we 
shall  integrate  over  the  whole  half 
space  to  compute  gains.  The  losses 
are  computed  from  beyond  a  certain 
radius.  Closer  to  the  sphere  we  must 


use  the  near  zone  approximation  to  be  discussed  below. 

Before  we  compute  the  energy  flux  we 
first  carry  out  the  calculation  mention( 
above.  In  order  that  a  particle  collid¬ 
ing  at  a  point  (r,  6}  hit  the  sphere, 
its  final  velocity  vector  must  be  in  the 
cone  shown.  Of  the  incoming  stream 
only  those  which  instantaneously  lie 

in  the  continued  cone  will  be  effective. 
The  ratio  of  the  projected  area  of  the  base  of  this  cone  to  the  total 
possible  impact  area  is  the  required  probability. 

A  CQsS-  , 

''G  ‘  "  Pr  =  9.  Q53  f 

A 


L 

\  /  J 

f— — — - 

P  -  g~ 


^*5 

'V 


‘irk!- 

=  =  TS*- 


r> 


The  factor  2  arises  from  inverse  collisions. 
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To  obtain  the  energy  carried 


the  eacolmter* 


per  particle  refer  to  the  geometry  of 


For  equal  masses  ^ 

M  =  i(T-e) 

'  f  .  C.I/  ^  j- 


(13) 


cj|  =  NQu 

NQ  =-L 


Call  q  =  energy  flux. 


± 

ij 


(14) 


We  fi^st  calculate  the  expected  loss  from 
Nj^(r)  we  must  take  1  (9)  s  | 


the  cylindrical  region,  hn  using 

(0)  -  X.  -  We  use  s  E  since  we  regard 
coUisions  as  carrying  away  all  the  energy. 

^hss  -  ^  ^ 


lAf 

1 


yji\ 


n 


Off 

e 


d 


f 


o<  -  ^ 

— T  .  '=<h 


The 


energy  gain  is  obtained  from  equations  (10)  through  (14). 

37r  ^  J ^(9)  QoUs^ny 

%*ss  ^  J  ^ 


(15) 


(16) 


(17) 
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There  is  then  a  slight  energy  loss  from  the  far  zone.  The  radial  integral 
will  be  needed  later,  it  is 


- 


r 


e 


L 


n. 


4  s  -C-  j:i_ 


£U  C"*  J 


(18) 


B. 


Near  Zone 


In  the  zone  close  to  the  sphere,  the  effects  ci  the  geometry  cannot 
be  glossed  over.  The  results  for  the  far  zone  indicate  a  net  loss;  a 
comparable  reduction  in  heat  transfer  is  found  to  take  place  in  the  near 
zone.  We  shall  present  a  first  attempt  to  estimate  what  occurs  in  the 

near  zone.  The  shaded  portion  in 
the  figure  is  the  region  in  which  there 
can  be  an  energy  loss  due  to  collision 
with  particles  coming  from  the  surface. 
However,  particles  which  otherwise 
would  not  hit  the  sphere  have  some 
chance  if  they  are  in  the  unshaded 
region.  A  particle  impinging  at  the 
stagnation  point  transfers  all  of  its  energy,  whereas  for  C  =  r/Z,  the 
maximiom  energy  transfer  is  about  half  the  available  energy.  The  energy 
given  to  the  sphere  will  be  taken  to  have  the  approximate  distribution: 


Er  -  ffl- 


■ )  ( '^/ry 


(19) 


The  density  of  reflected  particles  near  the  surface  is  largely  due 
to  their  low  velocity  and  to  the  cosine  distribution.  A  density  function 
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wMch  r;M3iild  approximate  physical  conditions  and  is  exact  on  the  surface 


The  value  cos  a  occurs  because  we  are  using  the  flux  weighted  component 
of  the  veiocityi  the  probability'of  a  particle  leavin’g  at  the  direction  d  is 
cOs  o|  its  contribution  to  the  flux  is  reduced  by  cos  a  and  finally  the 
component  of  the  velocity  is  u^  X  cos  a. 

,  To  compute  the  expected  gain  we  integrate  over  the  volume  between 


the  sphere  and  a  concentric  shell. 


)-ir 

m 


%  =  Nn-^f^0- 


The  loss  is  found  by  integrating  over  the  volume  indicated  by  the 


shading  in  the  figure. 


K.  It 


RcscS' 


fr  r ^  ^  3Tr  # 

°  »  P  J.  J  e 


Sin  R/r 


/o^S 


’  ^  R 
“  L 


C I  -  L-ft] 
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(25) 


The  ratio  of  loss  to  gain  is  then 


so  that  one  also  has  the  possibility  of  a  decrease  in  heat  transfer  in  the 
region  close  to  the  surface. 

We  now  bring  our  results  together  to  estimate  the  effect  of  the 

i 

region  in  between  our  two  extremes. 

(26) 


Near  Zone: 


‘gain 


Far  Zone: 


The  Jahnke- JEande  tables  list  values  of  the  exponential  integral 
needed  in  Eq,  (27),  For  small  values  of  the  argument  this  function  can 
be  calculated  from: 

tU-x)  =  -y  4  ^  X  -  (x-^) 

-Y  =  O.SIXX 

The  net  of  loss  and  gain  for  the  near  and  far  zones  is: 


‘near 


^far  ^ 


Q 


o  Wj,  Id 


“[-v‘  - 


(28) 
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lililll  IINIIiW  ^iiiiiaiwn*  . . . . 


The  net  from  the  entire  space  is  found  by  adding  the 
evalaatisg  the  result  at  r  s  33^, 
then 


net  energy  given  to  the 


above  and 
sphere  is 


C29» 


«>-PP«=dn.ation,  U.at  have  been  n^ade,  *e 

'«  t  the  factor  (0.  2)  is  known  to  better  than  t  SO*. 

A  correction  to  free  mol.cnl.  re.ult.  ntay  be  expected  in  Ute 
regime  where 


_a 

L 


In  terms  of  Reynolds  number  Eq,  0o)  is 


(30) 


(31) 


One  obtains  Eq. 


(31)  from  simple  gas  kinetic  arguments 


hf  =  Ma.ch  § 

^  “  Free  stream  temperature 
T^s  Body  temperature 


U 


k 


vp—  wkinematic**  viscosity 
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C-.  Numerical  Example 

In  order  to  have  a  feel  for  these  quanOTies,  consider  the  case  of  a 
satellite  at  an  altitude  of  100  miles. 


Some  pertinent  numbers  are 

Pa  ^  & 


a.4eoe  R/je:  ^ 


at  100  miles 


P/p^  = 


-  /a 
^  \o 


thus 


t  *■ 

|\1  s  IL-  6  ^  fo  Mokcofts  # 


%  L. 


The  air  molecules  are  regarded  as  rigid  spheres  of  diameter  3.  72A  ,  so 

that  -Is-  z. 

Q  10  CM 

/  '  N  Q  *  *7'  2-  ^ CM,"^ 


’flfZ.  c-H*  ^  oc- 

If  the  body  has  a  temperature  of  3^0°K,  then 

11  -  M-^TVIz  ..  . 


-3  U 

X/D  ^  ^  ~  0^*i 


Ufc  =  i 


m' 


—  H*  ^  1 0  C  Mi ,  A 


Under  these  conditions  we  have  the  following  results  based  on 


%o  --  NI(T7^“-F  =  X  ^  4a_ 


Eq,  {261 


I  fCS3 

^  0.0/ 9r 

po 

The  above  are  for  a  region  extending  one  radius  from  the  sphere. 
The  ratio  of  loss  from  the  region  greater  than  one  radius  away  from  the 
sphere  is: 


q  -  .  00’ 
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At  .00  .ph„. 

.  *  to  2%  t«.uc,..„  ,„  tou.  h«.  r..=ived  b.«us.  of  collision,  of  the 

ga.  wiU.  ^boondln,  n>oUcul„.  .bout  o,  .e  incident  nania« 

are  scattered  by  collisions.  However  th^  • 

nowe^er.  the  loss  is  considerably  less. 

since  many  of  the  deflected  particles  still  hit  tu 

P  rucies  still  hit  the  sphere  and  many  tha^ 

would  not  have  hit  are  deflected  into  the  sphere. 

From  Equation  (30)  we  have 

Vi)  ,  .  pffc 


—  3*  /  X  {0 


a 


(32) 


where  for 


6  h  c 

b  =  62 


ft. 


^  '  ^'30^ 

h  >/  /^o  itm. 

^  2.2, IT  f-f. 


Kio 


8 


we  conclude  that  for  a  100  cm  satelUte  traveling  at  26,  000  ft/sec 

and  ramng  a  surface  temperature  300°K,  free  molecule  flow  theory  gives 

a  re»  sonable  result  for  the  total  evT.>e. 

the  total  energy  imparted  to  the  sphere  as  long  as 

the  a  itude  is  greater  than  J 15  ?  10  miles. 
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ADDENDUM 


The  density  of  particles  as  a  function  of  distance  away  from  the 
sphere  r  along  the  direction  of  flight  (8  —  0)  is  given  by  the  expression 


■n  (r.^=e)  :=  i-fC- 


fy,. 


c,£>s  S' costx.  slfiS"  d  ^ 


This  quantity  is  plotted  in  the  accompanying  figure.  For  comparison  we 
also  plot  4R/r  and  4/3  -f  R/r.  The  former  is  our  fit  for  the  near  zone. 
The  latter  is  asymptotically  correct  for  the  far  zone;  howeverj  we  have 
only  taken  the  leading  term  which  is  4/3.  Our  fit  is  indicated  in  the 
figure  as  a  bold  line  and  can  be  seen  to  be  slightly  larger  in  the  near 
zone  and  smaller  in  the  far  zone  than  the  correct  density  which  is 
given  by  the  dashed  line.  It  is  felt,  however,  that  the  explanation 
should  give  a  reasonably  close  answer  to  the  problem. 
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ADDENDUM  II 


In  this  section  we  find  an  expression  for  the  energy  return  to  the 
sphere  which  is  exact  on  the  surface.  This  modifies  our  work  in  the 
near  zone  but  leaves  the  far  zone  results  unchanged.  The  new  result 
indicates  a  lower  coefficient  than  given  in  Eq  (29).  0.12  instead  of 
0.  18. 


A  collision  on  the  surface  can  send  in  both  particles  which  transfers 
all  the  available  energy  or  just  one  particle  which  carries  most  of  the 
energy.  The  former  occurs  if  the  final  velocity  vectors  of  both  particles 
are  directed  toward  the  sphere.  This  criteria  is  expressed  i.  i  terms  of 


the  tangent  plane  to  the  sphere. 
Consider  a  sphere  with  a  tangent 
plane  at  point  P  and  a  coordinate 
system  at  that  point  parallel  to  the 
old  one.  Measuring  angles  from  the 
-Z  axis  at  this  point  we  see  that  for 
two  particles  to  hit  the  sphere  that 
the  velocity  vector  of  one  of  them 
must  lie  between  zero  and  0 


m 

/ 

The  limiting  angle,  8^,  is  in  general  a  function  of  p  and  (S  • 

The  relationship  between  these  quantities  can  be  found  by  using  the  fact 
that  the  normal  to  the  tangent  plane  is  perpendicular  to  every  line  in  it. 

With  respect  to  the  coordinate  system  at  the  center  of  the  sphere, 
the  inward  directed  normal  has  direction  cosines  (-sin^  ,  0,  cos  <3>)- 
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In  the  same  coordinate  system  a  line  in  the  plane  will  have  direction 
cosines  (sin  6  cos  ^  ,  sin  6  sin  ^  ,  cos  6).  Since  these  lines  are 
perpendicular,  we  have 


sin©  sin  6  cos  +  cos  6>  cos  6  =  0 

tan  e  = 

cos  ® 


One  can  now  compute  the  number  of  particles  and  energy  returned 
to  the  sphere.  In  the  following  we  make  use  of  the  fact  that  for  a  collision 
such  that  the  velocity  vector  of  one  of  the  particles  lies  within  0  and  8^, 
two  particles  hit  the  sphere,  whereas  for  the  region  between  6^  and 
ir/2,  only  one  particle  goes  into  the  surface. 


The  number  of  particles  returning  to  the  surface  is 


where 


*  C  sa  Qj 


Adding  and  subtracting  a  in  the  numerator,  we  have 


The  integral  on  the  right  hand  side  is  listed  in  Groebner  and 


Hofreiter 


[,  #  331 . 56  a.  Our  result  is  then 
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=  1  +  COS0 


(3) 


The  appropriate  energy  return  function  is  found  from  an  equation 


similar  to  JEq.  (2) 

-  -L-  E 


'^C  i 


-  B 


TT 


^^ccsS-s,A^iS-  J  s.aJ- 

V2_ 


(4) 


The  integral  is  evaluated  in  the  same  manner  as  the  one  in  Eq.  (2). 


4  f'  ^ 


Si 


r 


IT/-. 


J 


AjL 


(cl'"  - 


The  first  integral  is  the  one  we  have  already  used;  the  second  integral 
is  proportional  to  the  derivative  of  the  first.  It  can  then  be  verified  that 
'^Itl 


Cos 


5,*A-  p') 


=  — f  i  -  CC5S  G  -  stP^gco^pj 


(5) 


and  that  the  energy  return  is 


=  -^P+CCSG  -t- ^  Ccj  0) 


(6) 


If  we  use  this  function  instead  of  Eq.  (19)  in  the  body  of  this  report, 
then  the  net  energy  loss  in  the  near  zone  would  be 
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The  outer  region  result  is  unchanged. 


The  total  net  evaluated 


at  r/R  =  3,  is 


total  net  energy  loss  = 


(0.12) 


(8) 


It  is  felt  that  the  "true**  coefficient  lies  somewhere  between  this 
result  and  (0.  2)  from  Eq.  (29). 
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